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%J . Abstract 



We present a class of singularity free exact cosmological solutions of Einstein's 
equations describing a perfect fluid with heat flow. It is obtained as generalization 
of the Senovilla class [1] corresponding to incoherent radiation field. The spacetime 
is cylindrically symmetric and globally regular. 
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Senovilla [1] has opened up a new vista of singularity free cosmological models. 
He obtained a class of cyclindrically symmetric inhomogeneous exact solutions of 
Einstein's equations describing a perfect fluid with the equation of state p = 3p. 
The spacetime was globally regular, smooth and satisfied the energy and causality 
conditions. All the physical and curvature invariants were finite everywhere, and 
all causal curves were complete [2] . There was no singularity of any kind. 

The singularity structure of a general class of inhomogeneous perfect fluid solutions 
was analysed [3] and it was shown that existence or non-existence of singularity as 
well as its kind (spacelike, timelike, Weyl, Ricci, etc) depended upon the particular 
choice of certain parameters and metric functions. There exists a general class of 
inhomogeneous perfect fluid solutions without singularity. Recently we have also 
found [4] a new class of solutions for stiff fluid (p = p) that have vacuum spacetime 
as the matter free (p = 0) limit. It is very curious that spacetime is empty as well 
as without singularity. One wonders about the source of its curvature. 

How robust is the singularity free character of spacetime ? In particular what 
happens when dissipative effects of viscosity, heat flow etc. are incorporated ? In 
this note we wish to obtain a generalization of the Senovilla [1] model that includes 
heat flow. 

Dissipative effects in cosmology were considered in the context of large entropy per 
baryon and isotropy of microwave background radiation[5,6] and heat flow was in 
particular investigated by several authors [7-16]. In very early times, matter is 
supposed to be highly dense and hot, and hence consideration of heat flow for very 
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early evaluation of the Universe is quite appropriate. 

The energy momentum tensor incorporating heat flow with perfect fluid reads as 



Tik = (p + p)uiUk - PQik + '^Ui^iQk) (1) 

where UiU^ = l,qiU^ = 0, and qi is the heat flow vector. The expansion 9, shear 
scalar a and the acceleration vector for the fluid are deflned by 



O-ik = U{i;k) - U(iUk) - ^(gik - UiUk) ■ 



ill = Ui-kU^ 



For Einstein's equations 



Rik — -^Rgik = -STrTik (2) 
we obtain the following class of solutions for perfect fluid with radial heat flow; 



ds^ = A^{dt^ -dr^) - B^dz^ - D^d(f)'^ (3) 

where 

A = C^-\at) C(^+4)/(i-2'')(3ar) 

B = C\at) C-(^+^)/9(3ar) (4) 

D=^S{Zar) C-(^+^)/9(3ar) C^^-''\at). 
Here a and b are free parameters and we denote C{x) = coshx,S{x) = 

sinhx, T{x) = tanhx, Se{x) = sechx. The metric (3) is globally regular for entire 
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range of the coordinates; — oo <t, 2;<oo, O<r<oo,O<0< 27r. The spacetime 
is cylindrically symmetric without any singularity. 

The physical and kinematical parameters have the following explicit expressions: 



a 



3(1 - 2b) 



{b + 4) (26 + 17) (2 - b)Se^{3ar) 



+ (5 _ b){b + 1) + a\l - b^)T\at) 



(5.1) 



8npA'={ ^^Y^) {b + 4n2b+17)Se\3ar) + ^{b + l){b + 7) 



W{l-b'')T\at) 



(5.2) 



87TqrA = ^^^—^{b + 4){2-b){b+l)T{at)T{3ar) (5.3) 

9A = a(2 - b)T(at) (5.4) 

a^A"" = la^{l - 2bfT^{at) (5.5) 
3n 

^r = -Yi:^(6 + 4)T(3ar). (5.6) 



The phenomenological expression for the heat conduction is given by 



q,=i;{F,i+Fui){Sl-u^Uk) (6) 

where i/j is the thermal conductivity and F is the temperture. For the case under 
consideration, the above equation can be integrated if -0 is a function of t alone. 
The temperature distribution in the Universe is then given by 
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F = lo{t)C^'+''^/^^-^''\3ar) - -^(2 - b){b + l)T(at) C'^^+^^/^^-'^'^Sar) (7) 

OOTTip 

where Iq is an arbitrary function of t. 

The solution(3) reduces to the SenoviUa [1] perfect fluid model for b — —1, when heat 
flux is switched off (q-^ = ip = 0). From the expressions in equation (5) the physical 
requirements {p > p> 0, p + 3p > 0) restrict the range for 6 to —1 < 6 < 1/2. All 
the parameters remain finite and well behaved for all values of t and r. We now 
compute the Weyl scalars and they are as follows : 



0? 



(1 - 2b) 



+ {l-b){l-2b)a^T^{at) 



(39 + 56 - 46^) 



2 / \ 2(i)+4) 

V'o - ^^4 = ^a\b -b)ib + A)Tiat)TiSar)C'^'-2^C ^'-\,^^^ 



g (at) (3ar) 



a' .,2^ 



1-26 



(35 + 136 - W) 



+ (1 - Zb + 2b')a''T\at) - + ^ J'' a^r^(3ar) 

(1 2oj 

It is clear from the above expressions that the Weyl scalars remain regular and 
finite for the entire range of coordinates for —1 <b< |. 

From relations (5) it follows that as t — >^ — oo all the physical and kinematic 
parameters, as well as the Weyl curvatures tend to zero, however the metric does 
not go over to fiat spacetime (this also happens in [1]). Then as t increases, the 
fluid starts contracting , the density and pressure increase reaching the maximum at 
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t = 0. The heat flux and expansion 9 are negative for — oo < t < 0, zero at t = 
and positive for < t < oo. Similarly the shear is positive in r and directions and 
negative in z direction and it reverses its sense at t = 0. The radial acceleration 
is zero at r = and negative for r > 0. As t increases from zero, the fluid starts 
expanding and heat flows radially out. All physical and kinematical parameters 
tend to zero again as t — > oo and the spacetime tends to the same state as that of 
t —>■ oo. The universe therefore starts from nearly flat spacetime at t — > — oo, passes 
through the dense and hot state at t = (the free parameter a can be chosen as 
large as we please) and becoming nearly flat and cold again as t — > oo. 

We have thus shown that it is possible to incorporate heat flux in the Senovilla 
class [1] of models without disturbing their singularity free character. We have also 
subsequently generalised the general class of perfect fluid solution [3] to include heat 
flux [17]. 
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